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The diffraction of electromagnetic waves by a convex cylindrical surface is considered. 
Attention is confined primarily to the region near the light-shadow boundary. The complex- 
integral representation for the field is utilized to obtain a correction to the Kirchhoff theory. 
Numerical results are presented which illustrate the influence of surface curvature and 
polarization on the diffraction pattern. Good agreement with the experimental results of 
Bachynski and Neugebauer is obtained. The effect of finite conductivity is also considered. 

1. Introduction 

The Kirchhoff theory of diffraction lias been quite successful in the prediction of field 
strength when the transmission path is obstructed by obstacles such as mountainous ridges [1-8]. 2 
In its usual form, however, this theory disregards the difference between vertical and horizontal 
polarization. Furthermore, the electrical properties of the diffracting obstacle are not con- 
sidered since in effect, the obstacle is regarded as a "black screen" or absorbing knife edge. 
A brief mention of the various developments, both old and new, will be made in order to properly 
orient the reader. 

Schelleng, Burrows, and Ferrell [1] in a now classic paper, applied knife edge Kirchhoff 
theory to the computation of the fields behind mountainous obstacles and other pronounced 
terrain features. The influence of the terrain in the immediate vicinity of the terminals was 
accounted for by the introduction of ground-reflected rays which are also diffracted by the 
obstacles. 3 A similar approach was made by Dickson, Egli, Herbstreit, and Wickizer [5], who 
demonstrated that the presence of an obstacle can, in some cases, produce an increase of field 
strength beyond that expected in the absence of the obstacle. This phenomenon has been 
described rather aptly as "obstacle gain" and has been used to advantage in the siting of the 
terminal and repeater stations in microwave relay links. 

The Japanese workers [4,6] have also been very active in this field. Presumably, the 
mountainous features of their country have motivated much of this work. Matsuo [4] in 
particular, has applied the Kirchhoff theory to many geometrical forms including conical shaped 
obstacles. wSome fundamental theoretical work has been carried out by Furutsu [9] who has 
presented solutions for diffraction by very general forms of obstacles. 

S. O. Rice [11] has given a comprehensive solution for the diffraction of a plane wave by 
a perfectly conducting parabolic cylinder which was considered to be a more realistic model 
of a terrain obstacle. Particular attention was paid to the fields near the horizon plane on the 
boundary region between light and shadow. It was shown that to a first approximation, the 
diffraction pattern behind the obstacle was very similar to that for a knife edge. The main 
effect of the finite radius of curvature of the crest was to increase the field strength for vertical 
polarization and to decrease it for horizontal polarization. A similar observation had been 
made by Artmann [12] sometime earlier who employed a circular cylindrical model. 

Still another approach to this problem was adopted by Norton, Rice, and Vogler [13]. 
They represented the diffracting obstacle first by a sphere for the case of a very large radius of 
curvature and then by a knife edge for vanishing small radii of curvature. A curve of diffrac- 
tion loss as a function of radius of curvature was then drawn by interpolating between these 
two limiting cases. According to their results, the diminishing of the radii of curvature is ac- 

1 A preliminary version of this paper was presented at the International Conference on Wave Propagation, Liege, Belgium on October 8, 1958. 

2 Figures in brackets indicate the literature references at the end of this paper. 

3 A scheme which will suppress or modify these ground-reflected rays has been suggested by TI. E. Bussey (Reflected Ray Suppressions, 
Proc. IRE 38, 1453, 1950). 
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companied by enormous increases in the field strength even though the separation between 
terminals and the angular distance is unchanged. 

Refinements of Kirchhoff diffraction theory when applied to the diffraction by smooth 
cylindrical models were made by Neugebauer and Bachynski [14, 15]. The basic idea of their 
analysis is that the equivalent radiating aperture is illuminated not only by the incident rays, 
but from rays reflected from the illuminated side of the cylindrical crest. The subsequent 
radiation from this aperture plane, which extends from the apex of the crest upward to infinity, 
is again reflected from the shadow side of the crest. Neugebauer and Bachynski are quite 
willing to admit that this is not a rigorous procedure, but nevertheless, it leads to results which 
are in good agreement with experiment. 

It is the purpose of the present paper to present a method for the computation of diffracted 
fields by convex surfaces which is both rigorous and easy to apply. The mathematical basis 
of this method was described in a previous paper [16]. This was an extension of the funda- 
mental work of Van der Pol and Bremmer, and Fock. The reader is referred to the above 
quoted paper for an account of the early theoretical work. 

Unfortunately, the existing solutions for the diffraction by spheres and cylinders are 
poorly convergent in the line-of-sight region. Consequently, these formal solutions require 
a great deal of numerical work despite the fact that they are mathematically elegant. A sig- 
nificant break-through was made when it was suggested that the complex-integral representa- 
tion of the field was easier to evaluate by numerical means than attempting to sum the residues 
of the poles of the integrand. While this point was probably recognized by Van der Pol and 
Bremmer [17], in the 1930's, it was not until 1946 that Academician V. A. Fock actually demon- 
strated this in what is now regarded as a classic paper [18]. Fock's method was applied by the 
present authors to the computation of the radiation pattern of an antenna mounted on or near 
a convex lossy surface [16]. The present problem is really an extension of this work. 

2. Formulation 

The problem is to calculate the field at some point P exterior to the cylinder in terms of 
the field on its surface (see fig. la). Since the incident magnetic field H inc (p, <£) is taken to 
have only a z component, the resultant field H (p, <j>) also has only a z component. An imme- 
diate application of Green's theorem yields 




Figure 1. (a) The infinite cylindrical obstacle, (b) Waves 
emanating from the penambral regions, (c) Geometry for 
source at A. and observer at A 2 . 
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n 2 {p,4,)=H™( P ,<l>) +af* [i7 2 ( P ',4>') ^ Q( P ,<t>; P ' ,*1-Q(p,<t>;p\<t>')(^~^)'] <fr', (1) 
where the Green's function is given by 

Q( P ,<t>; P ',<t>')=-lm»(kp), (2) 

with 

P=[p 2 +(p') 2 -2pp' cos (0-</>')] 1/2 . 

The primed coordinates p', </>', with p'=a, refer to a variable point on the surface of the cylinder 
whereas p, <j>, of course, refer to the position of the observer. The above integral formula can 
be simplified somewhat when the Leontovich boundary condition can be used [16]. Then 

dH ' (p / 0/) l =iwZHM,4>), (3) 

Op Jp'=a 

where Z is the surface impedance of the cylinder and e is the permittivity of the surrounding 
space. 4 Invoking this approximation, the external field is expressed explicitly in terms of the 
surface field Z/ z (#,</>) which can be regarded as known. The evaluation of the integral, how- 
ever, is not readily carried out unless certain restrictions are made. For the moment it will be 
assumed that the observer is at a large distance from the cylinder. Consequently, 

7 1 

f -ir/A 
/^_ e -ikp e ik p ' cos (<t>-<t>') f^\ 

4^Jkp 

where only the first term in the asymptotic expansion for H^ (kp) has been retained. 
In a previous paper [16], it was shown that 

H,(«,*^ft±/*' [G "*' )+! "" ] 9 [l-^+ 2 ™(^) , "] 

+h. ± .-* K'^**] , [y+i+2™ (!)"•], <5> 

for the region O^j |</>'| ^ (7r/2)+A</>, where (A<£) 3 << 1. The function g (X) occurring in the 
above expression is given by [16] 

1 /*oo -iXt 

where w x (t) and w[ (t) are Airy integrals defined in the appendix, and 

q= — i l — \ Z/r/o, with rjo^UOw. 

The quantity m can be physically regarded as the number of times a wave has crept around the 
cylinder. For (ka) l/3 7r>>l, only the term m = need be retained. Furthermore, if &a>>l 
and |0| less and not near tt/2 the geometrical optical approximation for H z (a, </>') is applicable. 
This reads 

H z (a,<t>')^H (l + R v e- ikacos +'), (7) 



* Equation (3) is strictly valid only when (Z) <<1207r bat actually this is more stringent than necessary [16j. 
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where R v is a Fresnel reflection coefficient given by 



„ _ (Z/ifo) — COS0' , , 

Uv (Z/i^+cos*'' W 



3. Evaluation of the Field 

The secondary or scattered field can thus be written in the following form 

fli( P) *) = ^_^(^+^), (9a) 



where 



f* CO 

ft= I [cos {4>'-4>) -kaZI Vo ]e- iM ^' 2)+ 4>'- cos i *'-^g(X l )d<t>' 



(9b) 



f 2 =\ [cos (<t>'-<t>)-kaZ/T,o]e- ikal( " 2) -*'- co ' ( <>'-'t> n g(X 2 )d<t>'. (9c) 

»/ — co 

The phase in the integral i/^ has the form 

&(</>') ^^[(Tr^+^-COS (<£' -<£)], 

so that the stationary points are at <t>' = cl>!i=(4n — 1) 7r/2 + for ?? = 0, 1, 2, 3 ... . The 
phase function is expanded about <j>' n and is of the form 

V(<l>')=ka[2Trn+4>+ 1 AW -$:)*], 

neglecting terms containing higher power than the third. The integrals over <t>' are now in 
the form of Airy integrals and thus fa and fa can be expressed in the form 

*>=^p .-a— «[-(!)"**} a* 

where 

where Wi (t) and v (t) are Airy integrals defined in the appendix. The quantity t is a new 
variable of integration. 

The special cases corresponding to q=0 and 5=00 lead directly to Goriainov [19] who 
assumed at the outset that the cylinder was perfectly conducting. These special cases are 
discussed in the subsequent text. The total field scattered by the cylinder is thus expressed by 

H'=-Ho£^(ho+hd, (12) 



where 



^^55^, (is) 
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and 




Ai=+V2(yJ {e- ika *G[{ka!2y'^]+e ika +G[-{kal2y^<t>]+e- 




+e- ika(2 *-+ ) G[(ka/2) 1 *(2<ir-<l>)]+ ...]}, 


(Ha) 


where 





0(*)=0(*)+£— _ (14b) 

The preceding results will now be generalized to include the case of a spherical wave 
incidence. The source is now imagined to be a magnetic dipole, located at (p , t, 0), thus the 
radius of curvature of the incident wave at the origin is p . Equation (12) for the scattered 
field in the transverse plane z=0 is now simply generalized to spherical wave incidence by 
multiplying it by Vpo/(p+Po)- This is valid for any form of cylindrical obstacle subject to 
p and p both large compared to the maximum transverse dimension of the effective scattering 
region. 5 Thus 

H --^(gp(^r (Ao+Al) ' (15) 

where H Q is the value of the incident field reckoned at the center of the cylinder and h and hi 
are as defined above. 

The portion of the scattered field proportional to h corresponds to Kirch hoff diffraction 
since it can be readily predicted by physical optics. It is independent of the electrical properties 
of the obstacle (and thus independent of the surface impedance Z) and is proportional to the 
cross-sectional dimension of the obstacle for small deflection angles. The remaining part of 
the scattered field, proportional to h h can be regarded as a correction to Kirchhoff theory. 
When tea,—*™, the ratio hi/h asymptotically approaches zero. For finite but large values of 
ka and small values of 4> only the first two terms containing G[(ka/2) l/3 <j)] and G[— (&a/2) 1/3 </>] are 
significant. These can be interpreted as waves emanating from the penumbral regions on the 
cylinder as indicated in the sketch shown in figure lb. In this case, the distances p and p are 
both assumed very large compared to the transverse dimension of the obstacle (i.e., 2a). To 
relax this condition, it is necessary to return to the original integral for the field arid use the 
following far-field approximation for the Green's function: 

4 v*p 

where p is not assumed to be large compared to 2a but is still large compared to the wavelength 
illustrated in figure lc. Carrying through the analysis leads to 



with 



c - ikR 

H^m+Hl^-jf- [F(X,u)+F(X',u')] ! (17) 

F(Z,«)=^fV ( « 2 rf«-^ e-'° 2 , (18) 

-\Jir J a U 

where 

and 

X=(ka/2Y^d. 

Equation (18) is only valid when u is somewhat greater than unity. A more accurate form is 

i c) 2 C(X) 
obtained when G(X) is replaced by G(X)~— 2 yzv 

5 See J. R. Wait, Appl. Sci. Research B.4, 464 (1955) for example. 
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The function F{X' , u') is the same form as F(X, u), and u' and X' are defined by 

and 

The quantities $i and s 2 are the linear distances from the source A x and the receiver A 2 measured 
to the intersection of tangent planes of the surface of the cylinder as sketched in figure lc. 
Similar remarks apply to the primed quantities. 6 The quantity R is the linear distance separat- 
ing A x and A 2 . 

It should be noted that when A x and A 2 recede to infinity, the earlier (Fraunhofer) solution 
is retrieved as a special case. To effect this transition, it is convenient to use the asymptotic 
relation 

[uF(X, u)] ^ e ~ ia2 Zjz G(X) for large positive a and (19) 

u-\ — j= G(X) for large negative a. 

Up to this point the analysis has been developed explicitly for the incident magnetic 
field vector directed along the axis of the cylinder. The results can be transformed to the case 
for the incident electric field vector directed along the axis of the cylinder by making the 
following transformation 

flS"(p,*)-*#?"(p,*) 

HI (p,cf>)-^E s z ( P ,<t>) 

[--«©■" (DH«~< &rm 

For metallic surfaces, q^O for II parallel polarization and g^oo for E parallel polarization, 
and these limiting cases are referred to in the following section as vertical and horizontal polar- 
izations, respectively. 

4. Numerical Results 

Basic to the present work is the function G (X). In the region near X=0, the relevant 
integrals cannot apparently be evaluated in closed form. The integrands are quite well 
behaved however, and it is not too difficult to evaluate them by numerical means. In figures 
2a and 2b the real and imaginary parts of G (X) are plotted as a function of X for q=0 and 
q=z oo ? respectively. 

Before discussing the diffraction pattern, it is of interest to consider the behavior of the 
field along the horizon plane. Clearly, this corresponds to 0=0, X=0, and a=0. The func- 
tion F (a) then becomes 

F(0)=~-^), (20) 



where 



/2fc« 1 « 2 \* (2/fca)* , , 

u=( — f-^ ) ' _ ? as before. 
\Si+s 2 / 2 



Numerical integration yields the following values: for q=0, G (0) = — 0.295+^ 0.0811; and 



« 6 and 6' are the angles subtended by the tangent planes as indicated in figure lc. They are referred to sometimes as angular distances 
and expressed in milliradians. 
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Figure 2. (a) Real and imaginary parts of the integral G(X) as a function of X for q = 0. (6) fled! aw/ imaginary 

parts of the integral G(X) as a function of X /or q= °° . 
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for q =co , G (0) = 0.344 — i 0.0918. This indicates that for vertical polarization (corresponding 
to #=0) with a metallic convex surface, the field along the shadow boundary is increased 
above that for a knife edge. The converse is true for horizontal polarization (corresponding 
to q — co). The increase or reduction from % is approximately the real part of the function 
[-0 (0)/u]. Therefore for 2 =0, \F (0)| ^0.500+ (0.295/u) and for q=<* j \F (0)|^0.500- 
(0.344/m). It is of interest to compare this departure from % with Rice's prediction using a 
parabolic cylinder model. The agreement is good if we equate the curvature of the apex of his 
parabolic cylinder with our a. It is also noted that his assumed plane-wave incidence is 
equivalent to moving our source to infinity (i.e., di = <x>). Then on identifying Rice's h with 
ka/2 and his A with [(si + S2)/(2irfcsiS 2 )]i, we find according to Rice for g=0, \F(0)\ ^0.500 + 
(0.296/w) and for q = co , \F(0)\ ^0.500— (0.342/w), being in close agreement with our values 
quoted above. On the other hand, Artmann, using a circular cylindrical model, gives for 
2=0, \F(0) | ^0.500+ (0.165/t0 and for q=™ , \F(0)\ ^0.500- (0.33/w). The disagreement in 
the case for q=0 with our result should be noted. It is suggested that the asymptotic approxi- 
mation for the Hankel functions used by Artmann is not fully justified. 

Neugebauer and Bachyinski have also computed a value on the shadow boundary using 
a semiempirical modification to Kirchhoff theory. Again transferring to our notation they 
would give 7 

271 
for g=0, \F(0) | ^O^OO-f-^- 1 -' and 

^8? 
for q= oo , \F(0) | ^0.500-^^- 

it 

The agreement of these values with ours (and those of Rice) is reasonable in view of the entirely 
different theoretical approaches employed. It is probable that at least part of the discrepancies 
can be attributed to the nonvalidity of some of the geometrical-optical arguments used by 
Neugebauer and Bachynski. 

5. Diffraction Patterns 

The diffraction pattern behind the obstacle is specified by the behavior of the function 
F(X, u) as a or uX varies. The amplitude \F (X, u)\ itself and the resulting diffraction loss L d 
are plotted in figure 3a for various values of u with q— and °° . As u tends to infinity, the 
patterns approach the Kirchoff form described by the function F (a). Since the function 
Fo(a) occurs frequently in what follows, its real and imaginary parts are plotted in figure 3b. 

The diffraction pattern behind a metallic cylindrical model of a mountain crest has been 
measured by Neugebauer and Bachynski in the K-band frequency range. For purposes of 
comparison, they also duplicated the experiment using a thin metallic sheet which was essen- 
tially equivalent to a knife edge. Their results for ka=239 are shown in figure 4a and 4b for 
vertical polarization (q=0) and horizontal polarization (g=oo), respectively. The ordinates 
are relative power and the abscissa are distance d 2 measured from the crest to the observer. 
The quantity h is the vertical displacement of the receiver above or below the horizon plane. 
In both figures 4a and 4b the source which is a transmitting horn is at a fixed location, being a 
distance d x (=150X) from the crest. Values calculated from eq (18) are also shown on the 
curves. The agreement is reasonably good. 

Neugebauer and Bachynski also measured the fields behind the crest for various values 
of the curvature parameter ha. This is shown in figure 5 where the ordinate is relative power 
and the abscissa is the deflection angle 9. The corresponding calculated curves are also shown 
in figure 5, taking ^^lSOX and d 2 =H3\. The agreement is still fairly good and probably is 
within experimental error. These curves illustrate in a striking way the influence of curvature 
on the diffraction fields. In the case of horizontal polarization and deep in the shadow, the 



i Their quantity 'A' = — 8w'/2 2-i/s times the real part of our G (0). 
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Figure 3. (a) Diffraction pattern for a conducting 
convex surface, (b) Real and imaginary parts 
of Fresnel integral F (a) as a function of a. 
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Figure 4. (a) Power distribution behind a cylindrical obstacle, (b) Power distribution behind a cylindrical 

obstacle. 



finite radius of curvature greatly reduces the fields below the Kirchhoff or knife-edge value. 
This is in accord with the results of Norton, Rice, and Vogler [13]. 

6. Further Numerical Results 

To obtain results for finite conductivity, it is necessary to evaluate the relevant integrals 
for finite values of ([. For the case of /^-parallel or vertical polarization 
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Calculated 

Measured (Bachynski 
and Neugebauer) 



_1_ 



Figure 5. Effect of curvature on diffraction by a 
convex surface. 



9, DEFLECTION ANGLE (MILLIRADIANS) 

where Z is the surface impedance. Now if the surface is a good conductor 

subject to a>>€co. That is, displacement currents in the conductor are neglected. The 
consequence of this approximation was discussed in some detail in a previous paper devoted 
to this subject [16]. Following the earlier notation 



where 



q^\q\e- iT/i =2 l '*Ae- iir/A , 



is the "loss factor" for the surface. It is seen that for infinite conductivity, A is zero. On 
the other hand, for very poor conductivity and large values of ka, the "loss factor", A, may 
approach infinity. These two limiting cases correspond to q = and q = co , respectively, and 
were discussed in the previous section. More usually A takes some intermediate value. 
For the general case 

G(X) = G 1 (X)+G 2 (X) 



where 



G(X) = G 1 (X) + G S (X), 

Jo' 



6 r I (X) = — ' v'{t)-gv{t) 



G 2 (X)- 



p-iir/4 

■yjT 

.Z37T/4 



w[(t) — qw x (t) 



dt, 



O 



Z7T/12 (*<*> 
\ 



G*(X), 



m=-r- 

m'tt Jo 



e (-Xt/2){^-i 



W 1 



(t) — qe- 2 * i/3 v(t) 



] 



dt. 



(21) 

(22) 

(23) 
(24) 
(25) 



The above relations follow immediately from eq (11) and (14b). The real and imaginary parts 

A A 

of the integrals G X {X) and G S (X) were evaluated numerically using tabulated values of the 
Airy integral functions. For the present purpose only the real and imaginary parts of the 
function G(X) are required. 

As mentioned above, the field along the horizon plana is of special interest. Expressed as a 
ratio to the free-space field, this can be written 
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/ 7 (o)= 



1 0(o) 



The real and Imaginary parts of G (0) are plotted in figure 6 in the form of an Argand diagram. 
Various values of A are shown on the curve. The left-hand portion of the curve (i.e., smaller 
A values) corresponds to an increase over the knife edge value. On the other hand the right- 
hand portion (i.e., larger A) corresponds to a decrease below the knife-edge field. It is of 
interest to note that for an A value of about 0.4 the magnitude of the diffraction field is almost 
the same as that of the knife-edge field. Associated with these intermediate values of A, how- 
ever, are appreciable phase shifts which correspond to a lag relative to the phase of the knife- 
edge field. The end points of the curve corresponding to q = and q=oo were discussed in the 
earlier section. 



Figurk 6. Integral G(O) 
plotted in the complex plane 
for various values of the con- 
ductivity parameter A. 

Applicable to vertical polarization. 
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The real and imaginary parts of G (X) are shown plotted in figures 7a and 7b, respectively, 
as a function of X for various values of A. 

The diffraction pattern defined by \F (X, u)\ can thus be readily computed for any value 
of A. Examples of these are shown in figures 8a, 8b, and 8c for u equal to 2, 5, and 10, respec- 
tively. 

The preceding results indicate that the influence of finite conductivity can be marked for 
//-parallel or vertical polarization. For purposes of illustration, the displacement currents 
were neglected so that q has a phase angle of —45 deg. For poorer conductivities this restric- 
tion may be violated and the numerical results are not strictly applicable although they can 
be used for approximate estimates as long as the conduction currents are at least greater than 
the displacement currents, (i.e., o->e w). In this case, the equation for A should be replaced 

The influence of finite conductivity for ^-parallel or horizontal polarization is negligible 
for most cases since \q\ is always very large compared to unity. 



7. Concluding Remarks 

A theory for diffraction of electromagnetic waves by a smooth and convex lossy obstacle 
has been presented. The model employed for computation is a circular cylinder of finite 
conductivity. 
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Figure 7. (a) The real part of the function G(X) as a function of X for various values of the conductivity parameter 
A. (b) The imaginary part of the function G(X) as a function of X for various values of the conductivity pa- 
rameter A. 

Both (a) and (b), applicable to vertical polarization. 
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^icure 8. (a) Diffraction pattern for a lossy 
convex surface u = 2. (b) Diffraction pattern 
for a lossy convex surface u = 5. (c) Diffrac- 
tion pattern for a lossy convex surf (ice u = 10 

(a), (b), and (<•), applicable to vertical polarization. 



Probably the most important conclusion is that the influence of finite curvature of the 
diffracting metallic obstacle is to add a positive correction to the Kirchhoff computed field for 
a vertically polarized incident field, whereas the correction is negative for horizontally polarized 
incident field. The increase or decrease of amplitude from the Kirchhoff result is accompanied 
by a corresponding change of phase which may be important in certain applications. 

The influence of finite conductivity is appreciable for vertical polarization and the curves 
presented in the paper may be used to compute this effect. In general, the ohmic loss tends 
to reduce the resultant diffracted field. In the case of a horizontal polarization, the finite con- 
ductivity has a negligible effect and the curves presented in the present paper for q=<*> may be 
used directly. 

While the present results are formally restricted to a diffracting surface of constant curva- 
ture, it can be expected that the results can be used for generally convex surfaces if the curva- 
ture is slowly changing. For example, it may be noted that the diffraction pattern has the 
same approximate analytical form for a circular cylinder, a sphere and a parabolic cylinder. 
At least this is true if the radius of curvature is large compared to the wavelength and the 
deflection angles are small. Deep in the shadow zone, this similarity principle breaks down 
and the field depends on the rate of change of the curvature. This aspect of the problem has 
been discussed recently by Wu [24] who has generalized some of the quasi-geometrical concepts 
of Keller [25]. As Wu points out, however, for small deflection angles (i.e., reasonably 
near line-of-sight) the diffraction field is not critically dependent on the analytical form of the 
surface provided it is sufficiently smooth. It can be expected that this will also be true for 
lossy obstacles. 
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8. Appendixes 

8.1. Definition of Airy Integrals 

The Airy integrals used in the above analysis can be represented by the following contour 
integrals 

W ^=J^ j mei2 ,n e a d*=u(t)-iv(t) (26) 

^((0=4= f i2 ,se«ds=u[(f)-iv'(t) (27) 

y IT J °° e% 

w 2 (t)=^= f ' e«ds=u(t)+iv(t) (28) 



y 7T i 

V 71 



wj(*)=-y= I se a ds=u'(t)+iv'(t) (29) 



where a = st — s 3 /3. 

The functions w(£), i?(£), u'(t), and y'(Z) have been tabulated by Fock [20] for real values 
of t 

8.2. Comparison With Results for a Spherical Model 

The previous discussion has been based on a cylindrical model. It is of interest to com- 
pare these results with those of Fock [21] who employed a spherical model. To facilitate this 
discussion, it is desirable to review briefly Fock's work as it pertains to this particular problem. 
The notation will be changed to be consistent with the rest of the present paper. 

The source is represented as a vertical electric dipole located at r^a+hi with respect to 
a sphere at r=a. (r, 0, </> are the usual spherical coordinates with origin at the center of the 
sphere and the axis in the direction of the dipole.) The great circle distance between the 
source and the receiver is denoted by s which, of course, is equal to ad. The receiver is located 
at r=a-\-h 2 . In Fock's theory the following parameter is consistently employed 

M=(kal2Y>\ x=M(s/a)=Md, 

7j l =kh 1 /M and y 2 =kh 2 /M. 

The fields for this geometry can be derived from a single scalar function U which has the 
following form 

U=e- ikR /R>+U s , 

where R is the distance from the source to the observer. The first term is a spherical wave 
and represents the primary influence whereas the quantity U s is the secondary influence 
resulting from the presence of the spherical body. The rigorous analytical expression for U s 
is well known, being given first by G. N. Watson in a form of a very slowly converging series. 
The transformation of this series into a complex integral had also been accomplished by Watson 
who suggested that a new series could be evolved if the contour was deformed round a set of 
complex poles. This results in the now well-known residue series which has been studied 
extensively. Rather than deforming the contour Fock decided that the complex integral 
could be employed directly. He approximated the spherical Bessel functions by an Airy 
integral approximation. This is equivalent to Van der PoPs and Bremmer's "Hankel approxi- 
mation'' which involved Hankel functions of order one-third. This representation is valid 
when the argument and the order of the spherical Bessel functions are very large, but nearly 
the same magnitude. This means that the heights hi and h 2 are small compared to the radius 
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a which itself must be very large compared to the wavelength. Fock is then able to write 
U in the following form 

p -iks 

u= e _ v 

■y/sa sin 6 

where V is an attenuation function which is expressed conveniently as the sum of two parts, 
such that V=V 1 -\-V 2 , where V 1 does not depend on the electrical properties of the sphere and 
V 2 is dependent on these properties. The explicit representation for these integrals are 

V -(^y 2 gj^ e- ixt w l (t-y 2 )w 2 (t-y 1 )dt+ J^ e-^w^t-yMt-y^dt^ (30) 

and 

*-- (tT \jI. '-" "M^§> -.«-*>»«-** 

+.fc, '-'" £w-£m «>.«-*>».«-!*><«} < ;il > 

where w { (t), w 2 (t), and v(t) are Airy integrals defined in appendix 8.1 and q is as defined previ- 
ously. The contour d is the straight line segment in the / plane from oo exp (— i2w/3) to 0, 
and C 2 is the segment along the real axis from to oo . By some ingenious manipulations, the 
above integral representation V l is shown by Fock to correspond to Kirchhon" diffraction and 
consequently, V 2 can be regarded as the correction to the Kirchhoff theory. More explicitly 



and 



V x ** f ^* e-^uF (a) for a>0 (32) 

ViM-*"* r^ 2 - , ^ x ttF (-a)l for a<0 (33) 



where oo ^k(R—$), ?/=(?/i?/2) 1/ 7(V?/i+'V / 2/2) 1/2 and a = u(jc—\ / y l —\!y 2 ). The function involves 
the Fresnel integral which is defined by 

F (a) = = — e-" 2 dx. 

y* J a 



It is now possible to express V 2 in terms of the integral G{X) encountered in the previous sec 
tion. In fact, 



V^-u^TT^-^GiX), (35) 



where X=a/u=x—^y l —^/y 2 . Equations (33) and (35) given above for \\ and F 2 , respec- 
tively, are valid if x, y u and y 2 are all large compared to unity, but in such a way that X or 
x—yjyi—^yo is finite. It is also important to remember that these particular results are only 
valid for h x and h 2 small compared to a. This fact appears to have been overlooked in a recent 
paper by Kalinin [22]. In many physical situations one is interested in computing diffraction 
by a spherical surface when hi and h 2 are comparable or even large compared to a. It is quite 
easy, however, to relax this restriction if the proper asymptotic form of the spherical Hankel 
functions h v ^ 1/2 (p) is employed and p is either a+h h or a+h 2 in Fock's work. First it should 
be noted that 

K. m (p)=J 1 ^HT{o), (36) 
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where H {2) (p) is the Hankel function of the second kind of order v and argument p. The Debye 
second order representation for the Hankel function Hf ) (kb) is now employed and leads to the 
approximation 

fl v-l/2(p)^^_ v 2l f) 2y/4> ( 37 ) 

where £=(p 2 — v 2 ) 1 ' 2 — v cos -1 (V/p), being valid when (p 2 — *> 2 )>>p 3/2 . With this change, it is 
now possible to write the solution in the form 

p — ikR 

^= V F, (38) 

where 

FsFoW+fJr'" 2 , (39) 

F,=-fl<(.X)/«, (40) 



VSi+*2/ 2 



_ a _f 2ks l s 2 \' 2 (2/ka) 1 * 

and 

X=(ka/2y'*6. (42) 

is the angular distance as usually defined, 8 and Si and s 2 are the linear distances measured 
from the source and the observer to the intersection of their tangent planes as in the previous 
section. 

The research reported herein was sponsored in part by the U.S. Air Force Cambridge 
Research Center, Air Research and Development Command. 
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